PPPL-2128 (Aug. 1984) 
Phys. Fluids 28(1), 116-126 (Jan. 1985) 



Efficiency of Current Drive by Fast Waves 

Charles F. F. Karney and Nathaniel J. Fisch 
Plasma Physics Laboratory, Princeton University 
Princeton, NJ 08544 



Abstract 

The Rosenbluth form for the collision operator for a weakly relativistic plasma is derived. The for- 
malism adopted by Antonsen and Chu can then be used to calculate the efficiency of current drive by 
fast waves in a relativistic plasma. Accurate numerical results and analytic asymptotic limits for the 
efficiencies are given. 
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I. INTRODUCTION 



Currents may be efficiently generated in a plasma by the injection of rf waves whose phase velocities 
are several times the electron thermal speed. 1 The efficiency, defined as the ratio of current generated to 
power dissipated, is achieved in this instance because the rf-generated plateau decays at a rate given by 
the collision frequency for the fast electrons, which is relatively low. In the quest for higher efficiencies, 
current drive by waves which interact with relativistic electrons has also been considered. 2 Relativistic 
effects modify the scaling of the efficiency, placing an upper bound on the efficiency achievable by current 
drive by fast waves. In this paper, we do several things: we give a more complete analysis of this problem 
based on a formalism adopted by Antonsen and Chu. 3 Specifically, we find that the effect of finite 
electron temperature leads to an enhancement of the efficiency. In order to calculate this effect, we first 
give expressions for the most important terms in the electron-electron collision integral in the relativistic 
limit. These expressions are put in Rosenbluth form so as to be amenable to easy implementation on a 
computer. We imagine that the relativistic Rosenbluth potentials that we identify may be useful in other 
problems arising in very hot plamas. 

In order to put the present work in perspective, let us briefly review the chief tools used in the study 
of current drive. The early work used fairly crude analytical models. 1,4 These models were sufficient 
to obtain the scaling laws for the efficiency of current drive, but were unable to provide the coefficients 
with any accuracy. Therefore, the analytical treatment was supplemented by numerical solutions to the 
two-dimensional (in momentum space) Fokker-Planck equation, 5-8 from which accurate estimates of the 
efficiency could be found. The first accurate analytical treatment of current drive was based on a Langevin 
formulation of the electron motion. 2,9 This involved taking the electron temperature to be small, allowing 
energy scattering to be ignored. The moment hierarchy for the Langevin equations can then be closed, 
which allows an analytical solution to be obtained. This was followed by a more complete numerical 
study of the Fokker-Planck equation for current drive in which the problem was reduced to the numerical 
solution of a one-dimensional integro-differential equation with a source due to the rf. 10 In this work 
toroidal effects were also included. The results agreed with the Langevin analysis 9 in the limit of large 
phase velocities (as they should) and gave more accurate numerical data for phase velocities comparable 
to or smaller than the thermal velocity. More recently, Antonsen and Chu 3 and, independently, Taguchi, 11 
using methods first used in the study of beam-driven currents, 12,13 recognized that it is not necessary to 
solve the rf-driven Fokker-Planck equation in order to find the rf-induced current. Instead, they showed 
that the Green's function for the current is the Spitzer-Harm function 14 describing the perturbed electron 
distribution in the presence of an electric field. This reduces the problem to the determination of a single 
two-dimensional function, from which the current generated by any form of rf drive can be calculated by 
a simple integration. 

Up until now, the only reliable analytical results for current drive in a relativistic plasma are those 
obtained using the Langevin methods by Ref. 2. As we will show, these are only exact for T e <C m e c 2 
and p 2 3> m e T e (where p is the momentum of the resonant electrons). A more complete analytical or 
numerical treatment along the lines of that achieved in the nonrelativistic case was hampered by the lack 
of a convenient form for the relativistic collision operator. This is remedied to some extent by the results 
of this paper where we calculate the collision integrals for the first Legendre harmonic of the perturbed 
electron distribution neglecting electromagnetic effects on the binary interaction (in this approximation 
the collision integral reduces to the Landau form). Having done this, we are able to generalize the 
treatment of Antonsen and Chu 3 to the relativistic case. A number of useful results flow from this: we 
can numerically calculate to high precision the current-drive efficiencies in the relativistic regime. We 
can perform an asymptotic analysis of the Spitzer-Harm problem to obtain analytic approximations to 
the efficiencies. In addition, we give higher-order asymptotic corrections to the current-drive efficiencies 
in the nonrelativistic limit. Throughout this paper, toroidal effects are entirely ignored. Although these 
effects are important in the study of current drive by low-phase-velocity waves, they play little role in 
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current drive by fast waves. Incorporation of these effects, however, proceeds in exact analogy with the 
treatment for the nonrelativistic case. 3 

Relativistic effects on rf current drive have also been considered by Hizanidis and Bers. 15 They take 
moments of the kinetic equation. In order to close the resulting system of equations, they approximate the 
steady-state electron distribution by a delta function. This approximation is unjustified and, consequently, 
their results for the current-drive efficiency are incorrect. 

The plan of this paper is as follows: In Sec. II we show how the relativistic collision operator may be 
reduced to the Landau form. In this form, the collision operator is costly to evaluate numerically. So, in 
Sec. Ill we convert the collision integrals to a Rosenbluth form, which may be evaluated very efficiently. 
The formulation of Antonsen and Chu is generalized to the relativistic case in Sec. IV. The numerical 
results for the efficiencies are given in Sec. V and the asymptotic results in Sec. VI. Finally, in Sec. VII, 
we examine the asymptotic form of the efficiencies using the full relativistic collision operator. 



II. RELATIVISTIC COLLISION OPERATOR 



The collision operator for a relativistic plasma is given by Beliaev and Budker. 16 They give the 
collision operator as 



dfa(p) C ' 



dt 



£C(/ 0) / 6 ), (la) 



where the kernel U is given by 

II 7a7fc(l ~ fla ' fib) 2 t r ,2/2/! a n.1 \2 ill 

- 7 a 2 AA - 1?0A + 7 a 2 7f (1 - /3 a • b )W b + 0A)}- (2) 

Here a and b are species labels, q s is the charge of species s, logA a / b is the Coulomb logarithm, e 
is the dielectric constant, p is the momentum, v s = cf3 s = p/m s -f s is the velocity of species s, and 
7 S = (1 +p 2 /tiIjC 2 ) 1 ' 2 . The distributions are normalized so that 



/ 



fs(p) d 3 p = n s 



the number density. We are primarily interested in situations where fast electrons are colliding off a 
weakly relativistic background. In that case [3' b <C 1, and we can approximate U by its nonrelativistic 
form 

u 2 \ uu , 

U = = , u = v a -v b . (lc) 

u 6 

Since the original form for U was symmetric in the primed and unprimed variables, we could equally 
well have obtained Eq. (lc) under the assumption that (3 a <C 1. The relative difference between Eqs. (lc) 
and (2) is 0((3' b ). However, the error in the collision operator C(/ , /&) is smaller than this. This point 
is examined in more detail in Sec. VII. Equations (1) are precisely the collision operator given by Lan- 
dau. 17 Indeed an examination of his derivation shows that the mechanics of the collisions are treated 
relativistically; the interaction, however, is calculated nonrelativistically assuming a Coulomb potential. 
Use of Landau collision operator implies a neglect of the relativistic (i.e., electromagnetic) effects on the 
binary interaction. What we have shown here is that such an approximation is valid provided at least one 
of the colliding particles is nonrelativistic. 
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It is readily established that Eqs. (1) conserve number, momentum, and energy (£ s = m s c 2 'f s ), that an 
_ff-theorem applies, and that the equilibrium solution is a relativistic Maxwellian / s (p) oc cxp(— £' a /T), 
where £' s = (£ s — Vd ■ p)/-\/l — v^/c 2 is the energy in a frame moving at v<j, and T and Vd are 
independent of the species s. 

Throughout the rest of this paper we will restrict our attention to an electron-ion plasma. We assume 
the ions are stationary and infinitely massive (m,i — ► oo). This allows us to express the electron-ion 
collision operator in (p, /j,) space (where [i = p|| jp and || and _L are with respect to the magnetic field) as 

where 

_ n e q\ log A e / e 
Ane 2 ' 

= g t logA eA _ _(k 
9e logA e / e ~ g e ' 

and we have assumed neutrality <7 e n e + g^rii = 0. In Eq. (3) and henceforth we will omit the species 
labels from all electron quantities. 



III. GENERALIZATION OF THE ROSENBLUTH POTENTIALS 

For computational purposes, the Landau operator is not the most convenient form for the collision 
operator. If the plasma is azimuthally symmetric, a two-dimensional integration must be performed at 
each point in momentum space. If the number of grid points is N x N, this requires 0(N A ) calculations. 
This requirement is dramatically reduced in the nonrelativistic case by expressing the collision operator 
in terms of Rosenbluth potentials. 18 Unfortunately, although the Landau operator can be used without 
change to describe the collisions in a relativistic Coulomb plasma, the Rosenbluth form no longer applies. 
(The derivation of the Rosenbluth form from the Landau form requires, for instance, that (d/dp) ■ U = 
— (d/dp 1 ) ■ U, a relation that only holds nonrelativistically.) 

However, because the kernel of the collision integral Eq. (lc) has the same form as in the nonrela- 
tivistic case, it is possible to borrow some of the techniques of Ref. 18. We convert the p' integration in 
Eq. (lb) to v' space, substitute a particular Legendre component for /(p')> an d manipulate the resulting 
integrals into the form 

J \v-v'\PMh(v')d 3 v' 

or 

J iv-vr 1 ^(/iW)dV, 

which may be evaluated in the same way as Rosenbluth potentials 18 (P^ is a Legendre polynomial). 

Here we give the resulting expressions for collisions off a stationary Maxwellian background, i.e., 
C(f, fm), and for collisions of a Maxwellian off the first Legendre component of a background, i.e., 
C(/ m ,/z/i). In both cases only electron-electron collisions are considered. These terms are all that 
are required for the solution of the Spitzer-Harm problem (giving the Green's function for the rf current 
drive) and they suffice for an accurate numerical solution of the two-dimensional Fokker-Planck equation 
as described in Sec. V. 

Beginning with the case of collisions off a Maxwellian, let us start by assuming merely that the 
background is isotropic /(p) = fo{p). The three-dimensional integrals in Eq. (lb) then reduce to one- 
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dimensional integrals giving 



p 2 <9/i 



where 



47rr r f p ?/ 2 r°° 1 



F(p) = 
B(p) = 



4tiT 
in 

4tiT 
3n 



P /o(P ) ~ d P 



uo 



V 

rP q„2 _ ,/2 



p'f (p')2v/c 2 dp' 



(4a) 

(4b) 
(4c) 
(4d) 



Specializing to the case / = fm an d using the relation df m /dp = —(v/T)f m , we find that 

F(p) = (v/T)A(p) 

and the steady-state solution to C(f, f m ) = is that / is a relativistic Maxwellian 19 with temperature T 



fmip) 



where 



ATrm 2 cTK 2 (e- 1 ) 

£ = mc 2 j, 
T 



-exp(-£/T), 



(5) 



e 



(6 = 1 corresponds to an electron temperature of 511 keV), and K n is the nth-order modified Bessel 
function of the second kind. 

For later use we define here a thermal momentum 



Pt 



'mT, 



a mean-squared velocity 



d'p = -v,', 

m 



i-^e + fe 2 + o ( e 3 ), 



a thermal collision frequency 



rar nq A m\ogh 
Vt = ^pJ = 47regp? ' 
and a collision frequency normalized to the speed of light 

T nq 4 log A 
c m 2 c 3 47re§m 2 c 3 ' 

These frequencies differ by a factor of two from those used in earlier publications. 1,2,5,6,8,9 Specifically, 
we have v t — vq/2 and v c = v/2. This means that all our normalized efficiencies are a factor of 
two smaller than in these earlier papers. (We made this change because the normalized Fokker-Planck 
equation in the high-energy limit now has a simpler form. This convention is also used by other workers 
in this field.) 
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For p >• p t , the indefinite limits in the integrals in Eq. (4) can be replaced by oo, giving' 



.20 




(6b) 



(6a) 



Note that the frictional force F(p) reaches a constant value as p —* oo. This implies, for instance, that 
an electric field smaller than Tv 2 /qTc 2 cannot produce runaways. 21 On the other hand, the pitch-angle 
scattering frequency B(p)/p 2 continues to decay as p — > oo. As the energy of the electron increases, 
its effective mass increases. It is then more difficult to deflect the heavier particle. In this limit, pitch- 
angle scattering is negligible compared with frictional slowing down. This is to be contrasted with the 
nonrelativistic case where the pitch-angle scattering frequency and the frictional slowing-downrate decay 
as 1/p 2 and the two processes are of comparable importance. 

The implication for current drive is that the efficiency of parallel wave-induced fluxes, say by lower- 
hybrid waves, approaches a constant. This can be seen as follows: Nonrelativistically, the efficiency 
increases as p 2 . Relativistic electrons, however, slow down faster because they are heavier, and they 
also do not carry more current when pushed in the parallel direction. Each of these effects reduces the 
efficiency by 7 ~ p; hence the approach to a constant. 

The other term we shall need is C(f m , /i/i). This term is rather harder to compute. We define 
hip) = fm(p)xi(p) and write C(f m , nfmXi) = pfnJixi)- Again, we reduce (this time after much 
algebra) the integrals in Eq. (lb) to one-dimensional ones to give 



The term in square brackets in the last integral matches that in the first integral except for the interchange 
of the primed and unprimed variables. The simplification of Eq. (7) was achieved, in part, with the help 
of the symbolic manipulation program, MACSYMA. 22 

Equations (4) and (7) are now in a computationally convenient form. Their evaluation involves the 
determination of a number of indefinite integrals (the unprimed variables should be factored out of the 
integrals for this step), and the multiplication of these integrals by various functions of p. If the distribu- 
tion functions are known on a grid of N points, then the computational cost is just O(N). Furthermore, 
the calculation can be arranged so that nearly all the computations vectorize. 23 The general solution of 
the linearized electron-electron collision operator C(f, f m ) + C{f m , /) = is 



where a, b, and c are arbitrary constants. With a = c = and b = p |, this provides a useful check on 
Eqs. (4) and (7) and their computational realizations. 




(7) 



/ = (a + b • p + c£)f, 
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IV. FORMULATION 



We now turn to the calculation of the rf efficiency. There are three steps involved: the specification of 
the rf current-drive problem, the identification of the Spitzer-Harm function as the Green's function for 
the rf-driven current; and the solution of the Spitzer-Harm problem. 

We begin with the specification of the problem. This is just a standard application of the Chapman- 
Enskog procedure. 24 The most important assumption is that the collisional time scale is much shorter 
than the transport time scale (the time scale for heating the plasma by the rf). This places some restrictions 
on the rf drive. However, these are usually not severe ones in the case of fast-wave current drive because, 
even if the rf is strong, there are few resonant particles and, consequently, the heating rate is small. 

The effect of the rf is to induce an electron flux 

df 

S = -D • / (8) 
op 

in momentum space, where D is the quasilinear diffusion tensor. 25 In the Chapman-Enskog ordering 
this is taken to be of first order. The zeroth-order electron distribution is given by setting the collision 
term C(f,f) + C(f, /j) equal to zero. The general solution is a Maxwellian Eq. (5) with n and T 
arbitrary functions of time and position. For simplicity we ignore the spatial variations. Since the rf 
drive is particle conserving, we may take n to be a constant. A drifting Maxwellian does not solve the 
zeroth-order system since the ions are taken to be stationary. 

The first-order equation is given by substituting / = / TO (1 + ip) with tp ordered small to give 

= ^ • S + i^M/ m | logT, (9) 

where 

C(f) = C(f, f m ) + C(f m , f) + C(f, fi) (10) 
is the linearized collision operator, and {£) is the mean energy per particle 19 

(£) = \j £fm(p) d 3 p 

The last term in Eq. (9) represents the heating of the Maxwellian. The equation for the time evolution of 
T is given by the solubility condition for Eq. (9), which is obtained by taking its energy moment. Since 
the linearized collision operator is energy conserving (recall that we take the limit mj — > oo, so that there 
is no energy exchange between electron and ions), this gives 

n d 4L=P 
dt 

where P is the power dissipated per unit volume by the rf 



S-vrf 3 p. (11) 



[There is another solubility condition given by the density moment of Eq. (9). This is automatically 
satisfied by taking dn/dt = 0.] The solution to Eq. (9) is made unique by demanding that f m i\) have zero 
density and energy. 
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In the nonrelativistic limit, Eq. (9) is the equation solved numerically by Cordey et al. w However, 
since we are usually interested primarily in the current density generated by the rf 



J = q J W||/ m V>d 3 p, (12) 

and the efficiency of current generation defined by the ratio J/P, we usually do not need to know the full 
solution for tp. 

The method for determining the current without solving for tp was given by Hirshman 12 and by 
Taguchi 13 for neutral-beam-driven currents and was introduced into the study of rf-driven currents by 
Antonsen and Chu 3 and Taguchi. 11 The key is to define an "adjoint" problem 

C(f m x) = -QV\\f m - (13) 

Again f m X * s required to have zero density and energy. This is the Spitzer-Harm problem for the per- 
turbed electron distribution function due to an electric field E = Tp||. Using the self-adjoint property of 
the linearized collision operator J ipC{f m x) ^ 3 P = / xC(fmip) d 3 p, it is readily found that 

J = J S- A X (p)d 3 p. (14) 

In this equation \ plays the role of a Green's function for the rf-driven current. The ratio of Eqs. (14) and 
(11) gives the efficiency 



p 



J S-vd 3 p 



An important special case is when the rf excitation is localized. Then it is only necessary to know the 
position and direction of the excitation to determine the efficiency 

where all quantities are now evaluated at the position of the excitation. If we compare this method with 
the Langevin method of Fisch, 2 we see that x is the mean-integrated current due to a group of electrons 
released at p at t = 

poo 

X(P) =Q (v\\) dt. 



o 

The power of these results is that the calculation of J/ P does not require a solution of Eq. (9) for the rf 
distribution tjj. On the other hand, Eq. (13) must be solved for the Spitzer-Harm function x- This reduces 
to the solution of a one-dimensional integro-differential equation, which may be accurately computed. 
Furthermore, in the nonrelativistic limit, it has been tabulated. 14 This method also substantially reduces 
the parameter space to be investigated numerically. The solution of Eq. (13) depends on two parameters 
only, Z and 8. In contrast, the solution of Eq. (9) depends on various parameters specifying the nature 
of the rf excitation (for instance, the direction of S, the minimum and maximum phase velocities, etc.) as 
well as Z and 0. 

In order to determine therf current-drive efficiency using Eqs. (15) or (16), we must solve the Spitzer- 
Harm problem, Eq. (13). The solution \ consists of only the first Legendre harmonic, so we substitute 
X(P) = MXi(p) mto E q- (13) giving 



1 2 ,, .0X1 vA(p)d X i 2B(p) + TZ/v 

?dp P A{P) ^ ^ Xl + + qV = °' 
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where A(p) and B(p) are given by Eq. (4), the electron-ion term is given by Eq. (3), and I(xi) is given 
by Eq. (7). The fact that the solution of Eq. (13) consists of only a single Legendre component constitutes 
an additional advantage to this method of determining current-drive efficiencies. The solution of the full 
rf problem given in Eq. (9) consists, in general, of many Legendre components. Often some truncation is 
performed in computing these numerically. 

Equation (17) may be solved by approximate analytic methods either by expressing \ as a sum °f 
Sonine polynomials 24,26 or by formulating the equation as a variational problem. 12 These methods have 
the disadvantage that they generally fail to reproduce the correct asymptotic (large p) form for \- This 
failing does not affect the calculation of the electrical conductivity significantly since in that case \ is 
integrated with a weighting factor proportional to f m . However, it rules out such methods for the study 
of rf current drive, since the efficiency may depend on the local value of \- 

This leaves us either with asymptotic methods, which we apply in Sec. VI, or with numerical methods. 
Numerical solutions to Eq. (17) have been given in the nonrelativistic case in Refs. 14 and 27. Here we 
use a simpler method that avoids most of the problems with the application of boundary conditions. We 
cast Eq. (17) as a one-dimensional diffusion equation by setting the left-hand side to dxi/dt and solve 
this diffusion equation until a steady state is reached. (The initial conditions may be chosen arbitrarily.) 
The integration is carried out in the domain < p < p max and the boundary conditions x(0) = and 
x"(Pmax) = are imposed. The diffusion equation describes the physical problem of the evolution of 
the perturbed electron distribution in the presence of an electric field and is therefore guaranteed to give 
the correct solution of Eq. (13) without having to worry about spurious solutions that diverge at p = or 
p = oo. Since this is a one-dimensional diffusion equation, it may be solved by treating the differential 
operator fully implicitly (the time step may be taken to be large). The integral operator I(xi) is treated 
explicitly and is recomputed after every time step. In the calculations shown here, the momentum step 
size was taken to be pt/50, the time step was taken to be 1000/^i, and the process converged (i.e., the 
relative change in \i P er step was less than 1 part in 10 10 ) after about 50 steps. 

In the following sections we will also need the function G(p) = Xi(p)/p so that x(p) = P\\G(p). In 
terms of G, the gradient of \ is 

d 

gZXiP) = G(p)P\\ +P\\G p (p)f>, 

where G p (p) = dG{p)/dp. 

V. NUMERICAL RESULTS 

The solution for \ is given as a contour plot in Fig. 1 for Z = 1 and = and 0.01 . From these and 
a knowledge of S, the direction of the rf-induced current can be determined. In the nonrelativistic case, 
Fig. 1(a), x rises ever more steeply as p is increased, giving the favorable p 2 scaling for the current-drive 
efficiency. 9 On the other hand, in a hot plasma, Fig. 1(b), the slope reaches a constant (the contour levels 
are equally spaced), leading to a limit in the efficiency of the current drive. 2 

Figure 1 also shows that the contours become vertical for p\\ small. This indicates that pushing 
electrons with small p« in the perpendicular direction (as with cyclotron-damped waves) is not effective 
in generating current. Pushing electrons in the parallel direction is effective, especially for small pu , since 
the denominator in Eq. ( 1 6) can be small. In general, when the contours of constant energy (p = constant) 
cross contours of constant x, the efficiency can be very large. 

Turning now to the numerical results for the efficiency, we begin with the case of a localized spectrum, 
Eq. (16). Although this situation may not be realized in practice, it is important because it can help us to 
determine the best current-drive schemes by showing where in velocity space to induce the flux. There 
are two major classes of fast waves that have been considered for current drive, namely Landau-damped 
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waves (e.g., lower-hybrid waves) for which S = p|| and cyclotron-damped waves for which S = pj_. 
Taking the limit p± — » 0, we have 

J _ G(p)+ P G p (p) 



P v 
J _ pG p (p) 



(18a) 
(18b) 



P v 

for Landau-damped and cyclotron-damped waves, respectively. The efficiencies are plotted in Fig. 2 for 
Z = 1 and 9 = 0, 0.02, 0.05, 0.1, and 0.2 (these correspond to T = 0, 10, 26, 51, and 102 keV). 
The curves for = in the two cases are given analytically from Eq. (24); they agree with the results 
of Ref. 2. This confirms the earlier analysis and shows that it is exact in the limit of T <C mc 2 and 

p 2 > mT. 

We next consider current drive by a narrow spectrum of Landau-damped waves. In this case, all par- 
ticles satisfying the Landau resonance condition u — k»vu =0 interact with the wave, and the quasilinear 
diffusion tensor is 

D oc S(to — fc||U||)p||p|| 
a 7<5(p|| - mw p 7)p||p||, 

where v p = uj/ku is the parallel wave phase velocity. Assuming that the electron distribution is weakly 
perturbed, we can take / = f m in Eq. (8) to give 

S oc 7«||/ m <5(p|| - 77w p 7)p||. 

When we substitute this expression into Eq. (15), we obtain 



J 



r ( G{p)+ ^}^ GM \ lfm{p)pdp 

Jpo P 



(19) 



P v f°° 

/ lfm{p)pdp 

where po = mv p / {l — v^/c 2 ) 1 / 2 is the minimum resonant momentum. This efficiency is plotted in Fig. 3. 
In the limit v p — > 0, the efficiency becomes large. This demonstrates that current may be efficiently driven 
by low phase velocity waves as was proposed by Wort. 4 

A similar analysis can be performed for a narrow spectrum of cyclotron-damped waves. The situ- 
ation is more complicated here because the electron cyclotron frequency depends relativistically on the 
momentum 28 and because relativistic effects distort the diffusion paths. 8 In addition, the variation of the 
diffusion coefficient with p± depends on the harmonic number. This means that the efficiency depends on 
three wave parameters u/kn, fl/kn (ft is the rest-mass cyclotron frequency), and the harmonic number. 
We therefore will only treat this case in the nonrelativistic limit. 

In the nonrelativistic limit (0 — > 0, p/mc — > 0), the efficiencies for both kinds of waves have been 
calculated by Cordey et a/. 10 and Taguchi. 11 They considered a narrow spectrum of Landau-damped 
waves for which the efficiency is given by the nonrelativistic limit of Eq. (19) and a narrow spectrum of 
cyclotron-damped waves for which the diffusion coefficient is 

Dccvj^ - u p )p±p±, 

where I is the harmonic number and v p = (u) — iO) / fc|| . Assuming that / = f m in Eq. (8), the efficiency 
for cyclotron-damped waves is 



P 



P CO 

/ (p 2 -p 2 o) l frn(p)G p (p)dp 
/ (P 2 -Pl) l frn{p)pdp 



= m\ , (20) 
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where po = mv p . (Here we consider only the fundamental cyclotron resonance I = 1.) In Fig. 4, we 
plot these efficiencies normalized to the thermal quantities together with the asymptotic results, Eqs. (31) 
and (32a). For mv p ^> pt, the efficiencies scale as v 2 as predicted by Fisch and Boozer. 9 The l/v p 
scaling seen in the Landau-damping case for mv p -C pt is obtained by taking the limit v p — > in 
Eq. (19) to give 

■]__}_ I D (P±)fm(p±)G(p±)p± dp ± ^ 
P v p J D(p ± )f m (p ± )p ± dp ± 

Here we have included an arbitrary dependence of Donpi. In Ref . 6, three different types of low-phase- 
velocity current drive were identified, namely by Landau damping, transit-time magnetic pumping, and 
Alfven waves. These methods differ in the forms for D(p±) 

!1 (Landau damping), 

(p± /pt) 4 (transit-time magnetic pumping), 

[2 - (p±/pt) 2 } 2 (Alfven waves). 

The case plotted in Fig. 4 is the first one (Landau damping). Evaluating the integrals in these cases gives 




where the coefficients C are given in Table I. The coefficients for Z = 1 should be compared with the 
(less exact) results of Ref. 6 obtained by a numerical solution of the two-dimensional Fokker-Planck 
equation where the constants of proportionality are given as 4, 6.5, and 6.5, respectively. The coefficient 
Cl has been determined analytically by Cordey et al. 10 to be 



The dependence on Z indicates that the current is unaffected by electron-electron collisions. This result 
may be derived by taking the momentum moment of Eq. (13). The electron-electron collision term then 
drops out (from momentum conservation) and the electron-ion term is proportional to the numerator in 
Eq. (21). 

The last numerical example is one in which we relax the condition that / = f rn in Eq. (8). This 
allows us to find the flux S that develops in the presence of high rf power. In order to determine S, we 
numerically solve the two-dimensional Fokker-Planck equation 

g-CL_(fl + £■□.££, (22, 

until a steady state is reached. The numerical collision operator is defined as 

c num (/) = c(f, u + c(f m , m/i) + c(f, h), 

where fifi is the first Legendre harmonic of /. The electron-ion term C(f, /,) is calculated using Eq. (3). 

In order to justify our handling of the electron-electron collisions, let us consider the linearized 
electron-electron operator C(f, f m ) + C(f m , /)■ The first term describes the relaxation of the tail par- 
ticles on the bulk and the second describes the concomitant heating of the bulk. The linearization is 
justified even with strong rf, as long as /(p) s=s fm(p) for £ ~ T. The linearized electron-electron 
operator conserves energy, and if this were used in Eq. (22), there would be nothing to balance the power 
input by the rf (there is no transfer of energy to the ions in the limit m, — > oo), and so a steady-state 
solution to Eq. (22) would not be possible. In Eq. (9), this is handled by allowing the temperature of the 
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Maxwellian to increase slowly with time. In the numerical code, we adopt a different approach, namely 
to modify the collision operator so that energy is lost in an innocuous way. The term responsible for the 
bulk heating is the second term C(f m , /). Let us write / in this term as a Legendre harmonic expansion 



oo 



Of the terms in this series, only one, the k = term, contributes to the bulk heating. (The energy 
moments of the other terms vanish.) Thus in order to lose energy we drop the term C(/ m ,/o). Of 
the remaining terms in the series, only the first, the m = 1 term, is of importance — it is responsible 
for ensuring conservation of momentum. Thus we retain only this term and approximate C(/ m , /) by 
C(f m , m/i) t0 give the collision operator C nU m- 

The collision operator C nU m has the following properties: energy is not conserved (thus allowing a 
steady state to be reached); momentum is conserved; and quantities such as the Spitzer-Harm conductiv- 
ity, which are given solely in terms of the first Legendre harmonic, are correctly given. To justify the way 
in which energy conservation is handled, we may check that the results are insensitive to the details of 
how this is done. One such check is given below where we compare the efficiency given by the numerical 
solution of Eq. (22), in which energy is lost, and that given by Eq. (15), where energy is conserved. 

We assume that the rf diffusion term in Eq. (22) is caused by high-power lower-hybrid waves whose 
phase velocities lie between v\ and vi- Thus we take 



where D(p) is chosen to be large enough to plateau /. [Here we choose D(p) = 10 i/ t pf /(l + p/Pt)-] 

Figure 5 shows the steady-state solution of Eq. (22) for Z = 1, 6 = 0.01 (T w 5keV), v\ = 0.4c = 
Apt/m, and V2 = 0.7c = 7pt/m (the parallel refractive index satisfies 1.43 < tin < 2.5). Using the 
numerical solution for /(p) and S(p), and the definitions (11) and (12), we obtain J = 3.74 x 10~ 4 qnc, 
P = 1.28 x 10- 3 TOncV, and J/P = 0.293 q/mcv c . 

This is to be compared with the result given by Eq. (14) with the numerically determined flux S(p), 
namely J = 3.77 x 10~ 4 <7nc and J/P = 0.296 q/mcv c . (The figure for P remains unchanged since this 
depends on S alone.) These two sets of figures are within 1% of each other. The excellent agreement 
illustrates two points: the approximations made in the numerical collision operator, namely, the neglect of 
the heating term C(f m , /o), has little effect on the results for the current-drive efficiencies (discretization 
effects are probably a greater source of error in these results); and the analytic result Eq. (15) can be used 
to obtain reliable figures for the efficiency for cases of strong rf. What is needed in the latter instance is 
an estimate for the rf flux S. This may be found from a numerical solution of a two-dimensional Fokker- 
Planck equation (as here) or from an approximate analytical solution. Some saving may be possible using 
this method in conjunction with a numerical code: since S reaches a steady state sooner than /, it may not 
be necessary to run the code so long in order to obtain a reasonably accurate estimate for the efficiency. 



We have seen that the efficiency of current drive may be expressed in terms of the solution of the 
Spitzer-Harm problem, Eq. (17). This equation may be approximately solved in the limit p 3> pt- We 
will begin with the relativistic case and later treat the nonrelativistic limit. We start by writing down the 
normalized form of Eq. (17) in the limit p » p t . We chose normalizations based upon q, m, c, and 
v c . Thus momenta are normalized to mc, xi to qc/v c , J/P to q/mcv c , etc. We use the same symbols 
to represent the normalized and unnormalized quantities. The coefficients A(p) and B(p) are given by 




for vi < P|| / (m-f) < V2, 
otherwise 



VI. ASYMPTOTIC ANALYSIS 
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Eqs. (6), suitably normalized. The integral term may be evaluated by replacing the indefinite limits in 
Eq. (7) by oo, giving when normalized 



M 8 3/2 ( H a (Q,Z) + H b (Q,Z) \ 

\ vp V 2 J 



where H a and H b are definite integrals of \i ( an d thus independent of momentum) that must be deter- 
mined numerically. In the limit 8^0, both H a and H b are finite. In normalized form with p 2 3> O, 
Eq. (17) reads 



ev 2 



d 2 Xi (v ! 3 2\d X i 
dp 2 \ 6 v"/ 3 p J dp 



vp z \ v z J \vp v z J 

The error in this equation is exponentially small. 

We now make a subsidiary expansion in small 9. In the limit 6^0, several terms in Eq. (23) drop 
out leaving 

i d X i i + z 

v z op vp z 

This may by solved with the boundary condition \i {p = 0) = to give 

This is the result derived using the Langevin equations by Fisch. 2 For integer values of Z, the integral 
may be expressed in terms of elementary functions. In particular for Z = 1 we have 

Xi = ( -^-7 )(vp- 2 log 7). 



Of particular interest is the efficiency for large p since this gives a limit to the efficiency of current drive 
by fast waves. If we let p> 1, the integral may be approximately evaluated to give 

Xi (1 + Z) \ogp. 

If we now take 8 to be finite, Eq. (23) cannot be easily solved. However, we may solve it in the limit 
p> 1. We achieve this by writing 

Xi^ap + (3\ogp (25) 

in analogy to the situation with = 0. Substituting this form of xi into Eq. (23) and balancing terms of 
equal order in p gives 

i + e^ 2 H b 



a = — 2 (26a) 

from the O(p ) terms and 

/3 = _(l + Z-39V?)a-8'/» ga (26b) 

from the 0(p~ 1 ) terms. When the rf excitation is localized, the current-drive efficiency is given by 
Eqs. (18) that, with xi given by Eq. (25), read 

J 13 

- ~ a + - (27a) 
P p 

■^P 1 -^ (27b) 
P p 
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for current drive by Landau-damped and cyclotron-damped waves, respectively. [The factor of 1 jv in 
Eqs. (18) is replaced by unity in the limit p — > oo.] Equation (27a) (with p replaced by p ) also applies 
for current drive by a narrow spectrum as given by Eq. (19). In the limit of p — > oo, the efficiency of 
cyclotron-damped current drive vanishes, while for current drive by Landau-damped waves J/P — > a. 
In order to determine this limiting efficiency, either Eq. (26a) may be evaluated using the numerically 
found value of H b (Q, Z) (see Table II) or else the equation may be expanded as a series in 6 to give for 
p — > oo 

p«i + |e + H b (o,z)e 3 / 2 . (28) 

H b (0, Z) is tabulated in Table III. 

We now turn to the solution ofEq.(17) in the nonrelativistic limit 9 — > 0. We shall still consider only 
the limit p p t . The limits here are nonuniform. Equation (23) was obtained by takings 3> Pt followed 
by — ► 0. Here we will take the limits in the opposite order. To do this, it is convenient to renormalize 
Eq. (17) using q, m, p t , and v t as the system of units. In this case, J/P is normalized to qjp t v t . Making 
this change of normalization and taking the limit 8 —> is equivalent to formally replacing O by unity 
and substituting v = p, 7 = 1, and V 2 = 1 in Eq. (23) to give 



p 



a 



dp 2 \J* ' p ) dp 



\(l + Z-^)xi + ^+P = 0, (29) 

5 J \ p z J p z 



where H(Z) = H a (0, Z) + H b (0, Z) (this is tabulated in Table III). For p > 1 (in this normalization this 
means p p t ), we may develop an asymptotic expression for \i as a series in powers of p. Balancing 
the terms in Eq. (29) from 0(p) (the leading order) to 0(p~ 4 ) gives 

p 4 9p 2 Hp 9 _ 2 

Xl ~ 5 + Z + (5 + Z)(3 + Z) + 2 + Z + (5 + Z)(3 + Z)(l + Z) + {P >' 

For localized excitation, Eq. (18) becomes 

J p = ^Tz + ( 5 + zK3 + z) + JTz + °^ (30a) 

J 3p 2 9 9p- 2 



P 5 + Z 1 (5 + Z)(3 + Z) (5 + Z)(3 + Z)(l + Z) 4) (3 ° b) 

for Landau-damped waves and cyclotron-damped waves, respectively. The leading order terms here 
(those proportional to p 2 ) are exactly those derived by Fisch and Boozer. 9 

In order to compute the efficiencies for current drive by a narrow spectrum of waves, it is necessary to 
carry out the integrations in Eqs. (19) and (20). The following asymptotic series is useful for this purpose: 

f-OO 

I exp(- \y 2 )y n+1 dy = exp(- \x 2 )[x n + nx n - 2 + n(n - 2).t"- 4 + •••]. 

J X 

For n even, the series terminates and is exact. The efficiency for current drive by a narrow spectrum of 
Landau-damped waves, Eq. (19) becomes 

J 4^ 6(6 + z) Hv- 1 

P = 5T^ + (5 + Z)(3 + Z) + 2T^ + ° ( ^ } ' (31) 
For a narrow spectrum of cyclotron-damped waves, Eq. (20) gives 
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for I = 1 and I = 2, respectively. The effect of the integrations is to change only the higher-order 
0(Vp) corrections to the efficiencies. The leading order terms are the same as for the localized excitation 
Eqs. (30). Equations (31) and (32a) are plotted in Fig. 4. These closely approximate the exact results for 

v p > 2v t 

VII. HIGH ENERGY LIMIT OF COLLISION OPERATOR 

In the previous section, we derived finite temperature corrections to the efficiency limit found in 
Ref. 2. However, the collision operator in the Landau form Eqs. (1) was derived by assuming that the 
background electrons are only weakly relativistic or that 9 < 1. We must check, therefore, that the finite 
9 corrections to the Landau operator do not affect the formula for the efficiency limit Eq. (28). 

The linearized collision operator Eq. (10) consists of three collision terms. Since in all practical cases 
the ions are nonrelativistic, the ion term C(f, fi) needs no correction. The term C(f m , /) contributes to 
the integral term I(xi) in Eq. (17). However, this resulted in a 0(Q 3 ' 2 ) contribution to efficiency limit 
Eq. (28), so that corrections to this term will be of still higher order. 

Therefore, we need only consider collisions off a Maxwellian electron background C(f, f m ). Fur- 
thermore, if 6 is small and if p » p t , we may take v' <C v in the full collision kernel Eq. (2) and 
approximate U by its Taylor expansion about v' = 0. By retaining terms up to second order in v', we 
obtain 

C (ff)- rd °* + V > f ■ 1 ^ .2 df 

< [T,Jm) T77T— < ,;, 77— • "7".' - 



2 8 Pj \ 3 k 8p k ^ 8v' k T 2 dv' m dv' m 4 8p k , 

where summation over repeated indices is implied and the superscript (0) is used to indicate that U and 
its derivatives are evaluated at v' = 0. Evaluating these coefficients gives 



(0) _ vSjk - VjVk 
jk ~ ,,3 

(0) 

jk _ *Vj 



U 



dv' k 



O U jk _ 2VjV k - ni\ U ik 



2dvL3vL 



(This calculation was carried out using MACSYMA. 22 ) If we compare these with the equivalent expres- 
sions using U from Eq. (lc), we find that only the term proportional to /3 4 is new. The high energy form 
of C(f, f m ) is given by Eq. (4a) with A(p) given by Eq. (6a), F(p) = (v/T)A(p), and 



l_IL M- — 



In other words, in the high-energy limit the electromagnetic correction only changes the pitch-angle 
scattering term. The new term has no effect on the asymptotic form for the efficiencies Eqs. (27) because 
it is smaller by (3 4 than another term in B, which had no effect. 

Connor and Hastie 21 also give an expression for collisions of high-energy particles off a fixed back- 
ground. The corrections to Eq. (6) that they obtain differ from ours. This is possibly because the back- 
ground distribution that they treat is only approximately Maxwellian. 

VIII. CONCLUSIONS 

We have considered the problem of current drive by fast waves in a relativistic plasma. Let us briefly 
review the approximations made. The major one is the reduction of the full collision operator to Landau 



15 



form. We show in Sec. II that this holds if the background temperature is small, T <C mc 2 . The 
corrections to the Landau operator in the high energy limit are derived in Sec. VII and are shown to 
be small. The second important approximation is the linearization of the electron-electron collision 
operator. This is accurate provided the rf strongly affects only electrons on the tail of the distribution. 
The subsequent analysis leading to the formula for the current-drive efficiency Eq. (15) is exact. In order 
to apply this formula, it is necessary to determine the rf-induced flux S from Eq. (8) and the Spitzer-Harm 
function \ from Eq. (13). 

We considered two methods for computing S: either to assume that / = f m in Eq. (8) (corresponding 
to linear damping) or to solve the two-dimensional Fokker-Planck equation, Eq. (22), numerically. The 
latter method may be necessary in the case of high rf powers and wide spectra. Note that the efficiency 
can be accurately calculated even if the S is known only approximately since Eq. (15), being an integral 
operator on S, is insensitive to small errors in S. Often, useful information can be extracted from Eq. (15) 
even with very limited information about S. If the rf spectrum is known, we can make some estimates 
(based on either numerical or approximate analytical solutions to the Fokker-Planck equation) of where 
in momentum space the flux is largest. We can then use Eq. (16) to give the efficiency. 

The Spitzer-Harm function \ can be determined by solving Eq. (17) numerically as in Sec. V. Since 
this equation is just a one-dimensional equation, there is little difficulty in obtaining arbitrarily accurate 
results in this way. This method can be regarded as exact. Alternatively, we found asymptotic forms for \ 
in Sec. VI. From this we can write down analytical expressions for the current-drive efficiency in various 
cases as given in Eqs. (27), (28), (30), (31), and (32). 

The primary application of this work is, of course, to maintain a steady-state toroidal current in a 
tokamak reactor. The viability of this scheme depends upon the amount of circulating power that is 
required. Thus, an accurate calculation of the current-drive efficiency, as well as an assessment of the 
best possible efficiency, are of crucial importance. 

When applying these results to the study of steady-state current drive in a tokamak, it is useful to 
convert the efficiency J/P to I/W where I = AJ is the total current, W = 2-kRAP is the total rf 
power, A is effective poloidal cross-sectional area, and R is the tokamak major radius. This gives 



The last two equalities give the conversion from the normalized efficiencies given in the figures and in 
Sec. VI to practical units. Figures 2, 3, and 4 contain scales in these units. 

The present work calculates the efficiency that can be expected from an arbitrary wave-induced flux. 
It is possible, therefore, to come to some very general conclusions about the best possible efficiency that 
can be obtained by driving currents with different waves. In particular, there is a limit, given by Eq. (28), 
to the efficiency of current drive with fast waves, such as lower-hybrid waves, that interact through a 
Landau resonance with relativistic electrons. These waves are, perhaps, the most likely candidate for 
current drive in a reactor. 

The present calculations also apply to other types of current drive, for example, relativistic electron 
beams. 29 Here, the efficiencies will be similar to those of Landau-damped waves. Care must be taken, 
however, in interpreting experiments on relativistic electron beams because the assumption of a steady 
state is generally inapplicable. 

The equations developed here apply to other forms of rf current drive. Some of these may be very 
efficient, more so than lower-hybrid wave-induced fluxes. For example, if low-phase-velocity waves 
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interact through a cyclotron resonance with fast electrons, the rf flux may be nearly parallel to the constant 
energy contours, at the same time that the collisionality of the resonant electrons is small. This gives very 
high efficiency, but, in practice, these waves are much more difficult to generate than are lower-hybrid 
waves. 

Settling the question of the highest attainable current-drive efficiency with fast waves should enable, 
we hope, tokamak reactor designers to assess the practicality of using waves to drive steady-state currents. 
There may, of course, be other effects that present difficulties, such as the accessibility of the waves or 
nonlinear effects. On the other hand, there may be effects, such as the bootstrap current, which could be 
helpful. 

Finally, we hope that the form that we derived here for the relativistic collision operator, which 
enabled us to solve for the relativistic Spitzer-Harm function, will be of use in other numerical problems 
dealing with collisions in hot plasmas. 
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Tables 

TABLE I. The coefficients for the efficiency for the three types of current drive by low frequency waves. 



z 


C L 


Cm 


C A 


1 


3.76 


8.49 


8.09 


2 


1.88 


5.17 


5.07 


5 


0.75 


2.55 


2.60 


10 


0.38 


1.42 


1.48 



TABLE II. Table of efficiencies J/ P for Landau-damped waves in the limit v p — > c. The efficiencies are 
normalized to q/mcv c . 



e 


Z = 1 


Z = 2 


Z = h 


Z = 10 


0.01 


1.04 


1.03 


1.03 


1.03 


0.02 


1.09 


1.07 


1.06 


1.06 


0.05 


1.25 


1.20 


1.17 


1.15 


0.1 


1.55 


1.44 


1.34 


1.30 


0.2 


2.19 


1.91 


1.70 


1.61 



TABLE III. The coefficients H a (0, Z) and H(Z). 



z 




H 


1 


13.69 


21.12 


2 


9.13 


13.51 


5 


4.94 


7.01 


10 


2.88 


4.01 
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FIG. 2. Efficiencies for localized excitation for (a) Landau-damped waves (parallel diffusion) Eq. (18a) 
and (b) cyclotron-damped waves (perpendicular diffusion) Eq. (18b). The different curves show the 
efficiencies for various values of the temperature 9 as indicated. In all cases Z = 1. The top scale gives 
the kinetic energy of the electrons. The right scale gives the efficiency for a plasma with n = 10 20 m~ 3 , 
log A = 15, and R = 1 m. 
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FIG. 3. Efficiencies for narrow Landau spectrum Eq. (19) as a function of the phase velocity v p . The 
curves correspond to the various values of 9. In all cases Z = 1. The top scale gives the parallel index 
of refraction n\\ — c/v p . The right scale gives the efficiency for the same conditions as in Fig. 2. 




FIG. 4. Efficiencies for narrow spectra of Landau-damped (L) waves and cyclotron-damped (C) waves 
(I = 1) for the nonrelativistic case 6^0 and Z = 1. Also shown as dashed lines are the asymptotic 
results Eqs. (31) and (32a). The right scale gives the efficiency for a plasma with n = 10 20 m~ 3 , T = 
lOkeV, log A = 15, and R = 1 m. 
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p„/mc 

FIG. 5. Contour plot of the steady-state distribution / obtained by numerically integrating Eq. (22). Here 
Z = 1, © = 0.01, v\ = 0.4c, v 2 = 0.7c. The resonant region is indicated. The contour levels are chosen 
so that for a Maxwellian they would be equally spaced with Ap = toc/30. 
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